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Abstract 

The modified versions of tlie bag model equation of state (EoS) are considered. They are constructed 
to satisfy the main quahtative features observed for the quark-gluon plasma EoS in the lattice QCD 
calculations. A quantitative comparison with the lattice results at high temperatures T are done in the 
SU(3) gluodynamics and in the full QCD with dynamical quarks. Our analysis advocates a negative 
value of the bag constant B. 
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I. INTRODUCTION 



A transition to the deconfined phase of quarks and gluons, the quark gluon plasma (QGP), is 
expected at high temperature and/or baryonic density (see, e.g., Refs. [jj and |2| and references 
therein). In the present study of the QGP equation of state (EoS) we consider the system 
with zero values of all conserved charges. This is approximately valid for the QGP created 
in nucleus-nucleus collisions at the BNL RHIC and even better for future experiments at the 
CERN LHC. Up to now the strongly interacting matter EoS could be only calculated from 
the first principles within the lattice QCD. These calculations are done for zero or very small 
values of the baryonic chemical potential. The QGP exists at high temperatures T > T^, 
where the critical temperature corresponds to the 1^* order phase transition in the pure 
SU(3) gluodynamics or to a smooth crossover in the full QCD. The main results for the QCD 
deconfined matter EoS can be illustrated by the Monte Carlo (MCHattice results (LR) for the 
energy density e(T) and pressure p(T) in the SU(3) gluodynamics [3|]. The quahtative features 
of the EoS at T > Tc can be summarized as follows. The pressure p{T) is very small at the 
critical temperature, p{T^/T^ « 1, and rapidly increases at T > . At high T the system 
reaches the ideal massless gas behavior p = £:/3, thus, s{T) = aT^ . However, the constant a 
which regulates the high temperature behavior is about 10% smaller than the Stefan-Boltzmann 
(SB) constant ass ■ Both e/T^ and Sp/T"^ approach the value a from below. The interaction 
measure {e — 3p)/T^, called also the trace anomaly, demonstrates a prominent maximum at 
T = 1.1 Tc . Note that these properties of the gluon plasma EoS are also valid in the full QCD 

The bag model (BM) [3] was invented to describe the mass spectrum of the hadron states. 
Soon after that it was suggested [8|] to interpret the bag constant B as the non-perturbative 
energy density term in the deconfined matter EoS. For several decades, the BM EoS has been 
used to describe the QGP (see, e.g., Ref. [Oj). In its simplest form, i.e. for non- interacting 
massless constituents and zero values of all conserved charges, the BM EoS reads: 

e{T) = asBT' + B, p{T) = ^ T' - B , (1) 

where e and p have a simple dependence on T modified by adding the bag constant B ( "vacuum 
pressure"). The SB constant in Eq. ((Tj) is asB = 7r^/30 (c/^ + 7dp/8), where ds and dp are 
the degeneracy factors for the massless bosons (gluons) and fermions (quarks and anti-quarks). 
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respectively. 

The main goal of the present paper is to study the modifications of the bag model EoS. We 
consider simple analytical parameterizations for the QGP EoS which include a linear and/or 
quadratic in T terms in the pressure function to satisfy the qualitative properties listed above. 
The quantitative comparison with the MC LR in the SU(3) gluodynamics and in the full 
QCD with dynamical quarks jsj will be done in Sections 11 and III, respectively. The Section IV 
summarizes the paper. 



II. GLUON PLASMA EQUATION OF STATE 

In this Section a quantitative comparison of the modified versions of the BM EoS is done 
with the MC LR jsl in the pure SU(3) gluodynamics. The e/T^ and ?>p/T'^ were obtained in 
Ref. jsl by extrapolation to an infinite continuous system. We take the MC values of these 
extrapolated functions e/T^ and 3p/T^ at the same T/T^ points where the interaction measure 
{e — 3p)/T'^ has been simulated on the finite lattice 32^ x 8 . To determine the parameters of 
different models discussed below we will minimize the sums of the square deviations at these 
T/Tc points for {e — eMc)/T^ and/or 3(p — pmc)/T^, where e, p are the model functions and 
£mc, Pmc are the MC LR. 

The recent lattice estimate for the pressure at very high temperatures T/T^ = 10^ is still 



about 3% below the SB limit [lO|. The lowest order perturbative calculations give {ctsb — ly) oc 
g'^iT) oc l/ln(T/A). The calculations within the perturbative re-summation scheme are 
comparable with the LR at T = (3 -^ 4) and suggest that the dominant effect of interactions 
is to turn massless quarks and gluons into weakly interacting quasiparticles. The quasiparticle 



approach of Ref. [12| (see also recent papers [IS] and references therein) treats the system of 
interacting gluons as a gas of non-interacting quasiparticles with gluon quantum numbers, but 
with mass m{T) which depends on T. The particle energy u and momentum k are assumed 
to be connected as u = [e + m2(T)]'/' . The energy density and pressure take then the 
following form jl2| : 



e{T) = HkHk ^ + B*{T) ^ eo{T,u) + B*{T) , (2) 

2tc^ Jq exp(u;/r) - 1 

P{T) = ^ / k'dk ^ - B*{T) = poiT^oo) - B*{T) , (3) 

OTT^ Jq u exp{uj/T) — 1 
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where the degeneracy factor d = 2{N^ — 1) equals 16 for the SU(3) gluo dynamics. The ther- 
modynamical relation, 

- P{T) = e{T) , (4) 

leads to the equation for the function B*{T), 

dB* _ _ Ao{T,u) dm 
dT m dT ' 

where Aq = Eq — Spo, and Eq, po defined by Eqs. f l2|51) are the ideal gas expressions for massive 
bosons. If the function m(T) is known one can calculate B*{T) from Eq. ([5]) up to an arbitrary 
integration constant B. For m = aT, where a is a temperature independent parameter, the 
function B*(T) derived from Eq. ([5]) equals to [l^ : 

B*{T) = B~ i Ao(T,a;). (6) 

One obtains the energy density ([2]) and the pressure ([3]), 

e{T) = aT' + B, p{T) = ^T' - B, (7) 

which has the form of the BM ([1]) with constant a equal to 

3d 



n=l 



— K2(na) + — Kiina) 
An 



K{a) asB ■ (8) 



The Ki and K2 in Eq. ([8]) are the modified Bessel functions. The constant cr in Eq. (JTj) includes 
the suppression factor For a — )■ 0, it follows k — )■ 1, and Eq. ([7j) coincides with Eq. ([1]). 

The function K,{a) decreases monotonously and goes to zero at a — )■ 00. 

The modified SB constant a < ctsb allows to fit the high temperature behavior of e{T) 
and p{T). In what follows the BM EoS ([7]) is considered with B and a being free model 
parameters. The LR [3] cover the temperature range (0.89 4.5)Tc. We consider the high 
temperature phase ("gluon plasma") at T > Tc, where is a point of the 1*^* order phase 
transition. To be precise, let us note that we use the MC LR for T > 1.02 to avoid the 
uncertainties at T = Tc where e{T) has a discontinuity in thermodynamical limit. The fit of 
the MC LR for 3p/T' gives a = 4.62 and B = 1.56 T^, and it is shown by the dashed line in 
Fig. [1^. One observes a correct behavior, 3p/T' = cr, at high T and an abrupt drop near the 
critical temperature, 3p{Tc)/T^ ^ . These features of p(T) are in a qualitative agreement 



with the LR. A quantitative agreement is however unsatisfactory. Moreover, the temperature 
dependence of e/T^ calculated by Eq. ([7]) with a = 4.62 and B = 1.56 appears to be in a 
qualitative contradiction with the MC LR (see the solid line in Fig. [1^). 
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FIG. 1: (Color online) The MC LR for the SU(3) gluodynamics. The e/T^ (squares) and 5p/T^ 
(triangles) are extrapolated to infinite continuous system [3] . The dotted vertical and horizontal lines 
correspond to T/T^ = 1 and to the Stefan-Boltzmann constant asB = Svr^/lS, respectively. The 
dashed hnes show Sp/T^ and sohd ones e/T'^ for the BM EoS 1^. a: The fit of 3p/T^ with the BM 
EoS (HI) gives a = 4.62 and B = 1.56 T^. b: The fit of e/T^ with the BM EoS dZD gives a = 4.73 
and 5 = - 2.37 . 

One can alternatively start from fitting the MC LR for the energy density function 6{T) with 
Eq. ([7]). Unexpectedly, one obtains a rather good agreement with MC LR for e/T^ admitting 
negative values of the bag constant B . The negative bag constant B = — 2.37 and a = 4.73 
needed in Eq. ([7]) to fit e/T"^ leads, however, to an incorrect behavior oi p/T^ (see the dashed 
line in Fig. Wp)- 



A modification of the BM EoS ([7]) was considered by Pisarski 15| : 



e{T) 



aT^ - CT^ + B 



p{T) 



a 



B . 



(9) 



A presence of the T^-terms in p{T) and e{T) has been further studied in recent papers 16|. 
For brevity, we will refer to Eq. (El) as the "C-bag model" (C-BM). The fit of the MC LR for 
p/T^ with the C-BM EoS (ED is presented in Fig. Hi. It gives, a = 4.92, 5 = - 0.13 
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and C = 1.8 T2 . One finds an agreement of the C-BM EoS (E]) with the LR for 3p/T^ . In 
particular, ^p{Tc)/T^ ~ . However, e{Tc)/T^ ~ 3, which exceeds the LR. 




FIG. 2: (Color online) The MC LR are the same as in Fig. 1. The dashed lines show 3p/T^ and 
sohd ones e/T^ for the C-BM EoS Q. a: The fit of Sp/T^ with the C-BM EoS ([9]) gives a = 4.92, 
B = - 0.13 T^, and C = 1.8 T^. b: The fit of e/T^ with the C-BM EoS & gives a = 4.69, 
B = - 2.64 T^, and C = - 0.28 . 

Trying to improve the quantitative agreement with the LR, one may start from fitting the 
e/T^ with Eq. ([9]). One observes indeed a better agreement for e/T^ with the parameters 
B = — 2.64 and C = — 0.28 T^, being very different from those found in the fitting of 
3p/T^ . These new values of B and C lead, however, to a qualitative disagreement of Eq. ([9]) 
with 3p/T^ LR, as shown by the dashed line in Fig. [2)d. 

Comparing the BM EoS (CD and C-BM EoS (El) with the MC LR we have faced the serious 
challenge. Very different values of model parameters, B for Eq. (^^, or B and C for Eq. (Q, have 
been found depending on whether we start from fitting 3p/T^ or from e/T^ . By admitting 
negative values of the bag constant B in Eq. ([7]) or Eq. IQ, one obtains a good fit of e/T^ 
in the whole temperature interval T > Tc, but finds a disagreement with LR for Sp/T"^, as 
seen from Fig. [lb and Fig. [2b. This finding looks contra-intuitive in view that the functions 
e{T) and p{T) are in the one-to-one correspondence to each other due to the thermodynamical 
consistency equation (jl]). In Fig. [3] we show the differences between the pressure functions 
p{T) calculated in the BM EoS ([7]) or in the C-BM ([9]), with parameters obtained from the fit 



6 



of e/T^, and the MC LR for pressure Pmc(T) . The difference of the pressures is divided by 
. From Fig. [3] one clearly observes a linear temperature dependence of {p — pmc)/T^ ■ The 





FIG. 3: (Color online) A difference of the model pressure p and the MC LR pmc divided by 
T^. a: Pressure p is given by the BM EoS dH) with a = 4.73 and B = - 2.37 . The Sohd 
line presents the linear function 3.94 T/Tc- b: Pressure p is given by the C-BM EoS ([9]) with 
a = 4.69, B = - 2.64 T^, and C = - 0.28 . The sohd line presents the linear function 4.48 T/Tc . 



thermodynamical relation (jl]) does connect the functions s{T) and p{T) . This connection is, 
however, not symmetric in the two directions. If the function p{T) is known, one finds s{T) 
from Eq. in a unique way. However, if the function e{T) is known, Eq. is the T^* order 
differential equation for the function p{T) . The general solution of this equation involves an 
arbitrary integration constant. This results in a linear in temperature term in the function 
p{T) . Thus, for e{T) in the form of Eq. ([7]), a general solution of Eq. (j4]) for p{T) can be 
written as follows: 

e{T) = (tT* + B , p{T) = -T^-B-AT. (10) 

3 

The term —AT with an arbitrary constant A corresponds to a gene ral solution of the homo- 



geneous equation Tdp/dT — p = as was noticed in Refs. [17|, ll8|]. For brevity we call the 
EoS ([TO]) the "A-bag model" (A-BM). The A-BM EoS 1^, in contrast to the BM EoS d?]) and 
C-BM EoS (Q, gives essentially the same values of the model parameters a, B, and A either 
one starts from fitting Sp/T"^ or from e/T^ . Figure H] demonstrates a good agreement of the 
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FIG. 4: (Color online) The symbols show the MC LR in the SU(3) gluo dynamics [3]. The lines 
correspond to the A-BM EoS ^ with a = 4.73, A = 3.94 T^, B = - 2.37 T^. The dotted horizontal 
lines in figures c and d correspond to p/e = 1/3 and = 1/3, respectively. a: The squares are 
e/r^, triangles 3p/r^, and circles 3,s/(4T'^). b: The interaction measure {e-3p)/T'^. c: The ratio 
p/e. d: The speed of sound squared c^. 

A-BM ([TO]) with the MC LR for the thermodynamical functions e/T^, 5p/T^, and 3s /{AT^) 
(where s = {e +p)/T is the entropy density), interaction measure, {e — 3p)/T^, the ratio p/e, 
and speed of sound squared, = dp/de. For the A-BM EoS (fTOj) . {e — 3p)/T^ does not depend 
on the parameter a whereas the entropy density s{T) does not depend on the bag parameter 
B. 
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Let us consider the EoS which includes both C and A T terms, 



e{T) 



CT'^ + B 



p{T) 



a 



- - AT 



B 



referred to as the AC-BM. The standard BM EoS ([7j) corresponds to A = C = 0, whereas the 
C-BM EoS (ED and A-BM EoS ^ correspond to A = and C = in Eq. (Ill]), respectively. 
A comparison of the AC-BM EoS (HI]) with the MC LR for the e/T^ and 3p/T^ in SU(3) gluon 
plasma leads to C jT^ « 1 . We thus conclude that the AC-BM ( ITT]) for the gluon plasma 
is reduced to the A-BM EoS f|TO]) . Note that the lattice study of SU(Nc) gluodynamics with 
Nc = 3, 4, 5, 6, and 8 colors performed in Ref. 22(| reveals that Sp/T'^ and e/T^ divided by 
the corresponding SB limits follow essentially the same curves for different Nc . 

A physical origin of the linear in T term requires further studies. In this connection we 
remind the famous problem of Gribov copies 19| and his suggestion of the modified gluon 



dispersion relation, u{k) = \/k'^ ^ M'^/k? , where M is a QCD mass scale. It was shown in 
Ref. 20j that at T >> M this dispersion relation gives the SB limit e/T^ = Sp/T"^ = asB and 



power corrections of relative order 
cut-off K phenomenological model 



/T^ for p/T'^ and for e/T^ . This also resembles the 



2l| where u{k) = k 6{k — K) , i.e. low-momentum gluons 



are suppressed but high-momentum gluons are effectively free. 



III. QUARK GLUON PLASMA EQUATION OF STATE 

The LR for the realistic equation of state - the QCD with 2-1-1 flavors (light u-, d- and 
heavier s-quarks) have been presented by "HotQCD" [4] and "Wuppertal-Budapest" 
collaborations. We compare the modified BM EoS with the latest LR 6|. The continuum 
estimates of the LR results for p/T^, (e — 3p), and in the temperature range 100 MeV< 



T <1000 MeV are presented in Table 5 of Ref. {g]. At the highest available temperatures one 
observes the ideal gas behavior, p = e/S, but constant a = e/T^ is about 20% smaller than 
the SB constant cxsb — 15.63. This behavior is similar to the case of SU(3) gluodynamics. At 
small temperatures, T = 100 -^ 140 MeV, the LR are expected to be smoothly connected with 
the thermodynamical functions of the hadron-resonance gas. 

In contrast to the pure SU(3) gluodynamics with a 1** order phase transition between glue- 
balls and gluons, the transition from hadrons to quarks and gluons is a crossover. This smooth 
transition takes place in the narrow temperature range, T = 150 -f- 200 MeV where the energy 
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FIG. 5: (Color online) The symbols are the MC LR for the 2 + 1 QCD EoS [6]. The vertical dotted lines 
show T = Tj=152 MeV and horizontal lines show the SB constant asB — 15.63 in figure a, p/e = 1/3 
and = 1/3 in figures c and d, respectively. The lines show the C-BM ^ results at T > Tj=152 MeV. 
The model parameters are a = 13.01, C = 6.06 T^, and B = —2.34 Tf. a: The squares are s/T'^, 
triangles Sp/T'^, and circles 3s/{4T'^). b: The interaction measure {e - 3p)/T^. c: The ratio 
p/s. d: The speed of sound squared c^. 

density increases strongly. Several characteristic temperature points of the crossover transition 
are presented in Ref. |6|: T=145(5) MeV at the minimum value of cl{T), T = Tj=152(4) MeV 
at the inflection point of {e — 3p)/T^, T=159(5) MeV at the minimum value of p/e, T = 
Tmax=^^^{5) MeV at the maximum of {e — 3p)/T^. Non of these temperatures is the critical 
one, and the model fit of the QGP thermodynamical functions does not depend too much on the 
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choice of particular starting point in the range T = 150 ^ 200. However, the hadron-resonance 
gas expected at low temperatures gives a concave shape of the interaction measure {e — 3p)/T^, 



while the 



jR jo] show a convex shape near the maximum at T = Tmaa;=191 MeV. Thus we use 
the LR [6| above the inflection temperature T > Tj=152 MeV in our model analysis^. This 
temperature will be also used to present the model parameters. 

We start with the AC-BM to fit the Sp/T"^ and e/T^ LR for the high temperature QGP 
phase. The best fit corresponds to negligible values of the linear temperature term, A/Tj^ << 1. 
Thus, in contrast to our analysis of the LR in the pure SU(3) gluodynamics, the AC-BM f lTTj) 
is reduced to the C-BM EoS (|9]) for QGP LR. The found model parameters are equal to: 
a = 13.01, C = 6.06 T^, and B = -2.34 T^^. A comparison of the C-BM with LR ^\ at 
T > 152 MeV is shown in Fig. [5l It demonstrates a good agreement of the C-BM (j9]) with the 
LR for the QGP. In particular, the model leads to the maximum position Tmax — 189 MeV and 
the value of {e — 3p)/T^^^ = 4 which are very close to the LR. 

IV. SUMMARY 

We have considered the modifications of the bag model EoS. They are constructed to satisfy 
the qualitative features expected for the QGP EoS. We make also the quantitative comparisons 
with the MC lattice results for the SU(3) gluon plasma and for high temperature equation 
of state with 2+1 dynamical quarks p. Our modification of the bag model equation of state 
includes the following features: a suppression of the Stephan-Boltzmann constant; linear or 
quadratic in temperature term in the pressure function; a negative sign of the bag constant. 
These features are needed to describe the lattice data. The best fit of the LR for thermody- 
namical functions in SU(3) gluon plasma are found within the A-bag model (fTOj) . This model 
corresponds to: e = aT^ + B, p = — AT — B. A linear in T term in the pressure func- 

tion is admitted by the thermodynamical relation (jlj) between e{T) and p{T). The expression 
for the energy density looks formally the same as in the standard bag model (I7j). A principal 
difference from the standard bag model is a negative value of the bag constant B . 

The quantitative comparison with the MC lattice results for high temperature equation of 



^ The precise matching of the hadron-resonanc e g as and LR is beyond the scope of this paper. The discussion 
of a possible procedure can be found in Ref. |23j . 
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state with 2+1 dynamical quarks [6[ shows the best fit of the thermodynamical functions for 
the QGP within the C-bag model (Q: e = aT^ - CT^ + B,p = - CT^ - B. This model 

also requires i? < to fit the lattice data. Note that a negative value of B found for the gluon 
plasma and QGP does not contradict to the bag model hadron spectroscopy [7] which requires 
i? > at zero temperature. 
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